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Abstract. - This paper examines the effect of cooling on disentanglement forces in polymers 
and the implications for both single chain puUout and polymer dynamics. I derive the explicit 
dependence of the distribution of these forces on temperature, which is found to exhibit a rich 
behaviour. Most significantly, it is shown to be dominated by large fluctuations up to a certain 
temperature To that can be determined from molecular parameters. 

The effects of these fluctuations on chain friction are analysed and they are argued to 
undermine the traditional melt-based models that rely on a typical chain friction coefficient. 
A direct implication for first principles calculation of viscosity is discussed. This quantifies the 
limit of validity of such descriptions, such as Rouse dynamics and the Tube model, and pave 
the way to model polymer dynamics around the glass transition temperature. 



Polymer chain dynamics in dense melts are described quite well by Rouse dynamics and 
the Tube model Q that have become the accepted paradigm. The latter assumes that a 
chain moves in an effective tube formed by nearby chains and against friction produced by 
disentanglement from the rest of the network. This mean field type of description becomes 
increasingly inaccurate upon cooling and loses validity altogether around the glass transition. 
While there have been attempts to improve the model by considering more than one chain (e.g., 
through constraint release |^), it is worthwhile to explore extensions to the one chain picture 
before rushing to more complicated models. A detailed understanding of chain dynamics over 
a wide range of temperatures is also required to interpret recently emerging experiments of 
single chain puUout [Q. A model for the latter has been proposed recently in which the 
chain's mean field environment is replaced by a detailed consideration of the local dynamics 
with a focus on the distribution of forces needed to disentangle non-permanent entanglement 
points (EPs). An intriguing result of reference is that at low temperatures (below the 
glass transition) this distribution is dominated by large fluctuations. This result raises several 
questions: (1) while it is evident that the fluctuations must decay with increasing temperature, 
it is unknown how, and at what rate, this happens, nor is it known to what temperature do 
the fluctuations persist. (2) is it possible to quantify the effect of the large fluctuations on 
rheological properties such as viscosity and so determine the range of validity of melt-based 
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Fig. 1 - A non-permanent entanglement: the black circle is an anchor perpendicular to the page plane 
around which the disentangling chain slides. 



models? The latter issue is essential as a basis for modelling polymer dynamics around the 
glass transition temperature. 

These issues are the subject of this Letter. The probability density function (PDF) of 
the disentanglement forces and its temperature dependence are derived. The PDF is found 
to exhibit a rich behaviour: At high temperatures it is sharply peaked around very small 
forces, as expected, but upon cooling the disentanglement forces increase and a broad local 
maximum develops. Below some temperature a gap appears at the low forces end which 
gradually increases until it becomes complete at very low temperatures. When this happens 
the PDF diverges at a value, /o, which is determined by molecular data. An algebraic tail then 
appears at the large forces end of the PDF, giving rise to even larger fluctuations. I discuss the 
implications of these results for viscosity calculations from molecular information and derive 
the friction force on individual polymer chains. The connection between the broadening of the 
PDF and appearance of large fluctuations in chain friction is used to quantify how melt-based 
models break down in this regime. 

We start by deriving the disentanglement threshold force distribution using the model 
introduced in In this model a non-permanent EP is regarded as a local wrapping of one 
chain around an anchor chain which is fixed by the surrounding network. Disentanglement 
occurs when the wrapped chain slides around the anchor under a pulling force, as shown in 
fig. 0. The activation energy required to move m segments one position clockwise, Et, = mE, 
is determined by the molecular interactions and for anchors consisting of one chain m = 2 is 
expected (e.g., as shown in fig. Q). The thicker the anchor the higher the value of Ej,. Tugging 
at bead 3 with an increasing force / in an arbitrary direction exerts a force (/ cos 0) along the 
2-3 link and the disentanglement force is the value of / when sliding starts. For simplicity, 
it is assumed that once a chain starts slipping it continues to do so until disentanglement is 
complete. Once disentangled, the portion of the chain immediately after it starts to uncoil, 
building stress on the next EP. This process continues until the entire chain is unraveled. 
At very low temperature, slipping occurs when the work needed to move the chain by one 
segment a is equal to the activation energy Eb, f cos0 = E^/a. The value oi Eb/ks (ks being 
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Boltzmann constant) is expected to depend only weakly on molecular weight. In reference Q 
the force distribution was derived at low (essentially zero) temperature. Here we investigate 
the effect of both the distribution of 9 and the thermal fluctuations on the distribution of /. 

Since equilibrated polymer networks are isotropic down to the scale of a few monomers 
the angle is uniformly distributed but only values of 6 in the range < 6* < Omax = 7r/2 — e 
need be considered. This is because: (i) 6* < leads to unwrapping rather than sliding, 
and (it) values of too close to 7r/2 are irrelevant because the forces required to overcome 
the threshold /o in this regime are too large 1/6*) and would cause scission before sliding 
occurs. It follows that 



CO^Qmax — Eb/ {a f scission) = £ ^ 1 (1) 

where /scission is the force to rupture the chain and a is Kuhn's length. As the chain is 
pulled, the force on the EP steadily increases and disentanglement occurs when the energy 
input is larger than {Ef, — af cos 9). The probability of disentanglement to occur at all is 
therefore proportional to Q~l3(Eb-af cos e) ^]-^gj^ fcos6 < /o and identically unity otherwise. 
The probability of the former scenario is 



7(f) 



f 

1- / e 





and the probability of the latter one is 



fVif/for 1 



1 



-/3(iJi,-a/'cos(7(/)) 



df 
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In these and the following expressions (3 = l/fcsT, 7(/) = arccos(/o//), S{x) is Dirac's delta- 
function and H{x) is the step function which is 1 for x > and otherwise. The PDF that 
disentanglement occurs exactly when the pulling force reaches a value between / and f + df 
is therefore 



where 



Pif,T)df = C{Pi+P2)df = C 
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df 
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and the normalization constant is 



^l3af cos 9 
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-PEt 



d9 




When T ^ C — > 1/ 9max and the PDF converges to the first term, reproducing in this limit 
the result of reference . The second term gives the probability density of thermally assisted 
disentanglement and governs the range f < fo- Fig. |^ shows P{f,T) generated numerically 
at T = 0.05, 0.1, 0.15, 0.2, 0.4 and 0.9, where the temperature is measured in units of Eb/ksT. 

At T = 0.9 the PDF is narrowly peaked around / k, 0.05/o and as the temperature 
decreases it develops a maximum which moves towards higher values of /. Around T « 0.2 
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Fig. 2 - The PDF of the disentanglement forces P/(/) for T = 0.05(+), O.l(o), 0.15(*), 0.2(+), 0.4(|— |) 
and 0.9(x) pEb. Note the broadening and the development of a local maximum, which moves to 
larger forces with decreasing temperature. 



the PDF is smeared almost over the entire range < / < /o with a broad hump around 
/ ~ 0.5/o. The maximum sharpens again with further cooling while moving to higher values 
of /. Another notable feature is the appearance of a gap at small values of / as the temperature 
drops below T = 0.2, corresponding to an exponentially low probability of disentanglement 
by small forces. Around the same temperature a tail starts growing at values of / larger than 
/o. Finally, as T approaches zero, the gap is complete in the sense that no force below /o can 
cause disentanglement, the PDF is sharply peaked at /o and the tail assumes the algebraic 
form of the first term in eq. (^). Fig. ^ shows the dependence of the mean force, F = (///o), 
on temperature for = 100. In principle, {f(T)) should be computable directly from the 
PDF but the resulting expressions are cumbersome. Empirically, the data can be fitted very 
well empirically by 

(///o> - AT-e-y^ + B , (5) 

where A ^ 9±l, x = 0.57 ± 0.03, y = 11.2 ± 0.3 and B = 0.037 ± 0.002. It is important to 
note that the values of these parameters depend on molecular data through the dimensionless 
quantity e = E},/ {a f scission)- This has the intriguing implication that changing the chemistry 
to vary Ejj and f scission can have strong effects on the dynamics and rheology in the manner 
described by the results below. The very low-temperature behavior {T < 0.1) is dominated 
by large fluctuations of the form 

VWT)-^^+0{e-P^^) . (6) 

As the temperature increases the moments of / get smaller, the second term in eq. becomes 
dominant, and the sensitivity to e is washed out. Both the long tail at low temperatures 
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Fig. 3 - The mean disentanglement force as a function of temperature. The points are obtained from 
a numerical calculation and the line is a power-exponential fit, eq. Note that this curve is also 
proportional to the viscosity ^, eq. (W). 



and the broad distribution at intermediate temperatures suggest that large fluctuations are 
dominant below some temperature Tq. Note also that the smaller the value of e (the higher 
the value of /scission) the larger the fluctuations and therefore the higher is Tq, suggesting a 
way to manipulate Tq via the chemistry. 

To estimate the crossover temperature Tq consider the relative width of the PDF, F = 
{Spy^'^/if)- When F < 1 fluctuations are insigniflcant and the force is well characterised 
by a typical value, ftyp- When F is of order 1, fluctuations dominate and this is not a good 
approximation. Figure 4 shows that F increases steadily with temperature up to T w 0.4 
and then saturates to a plateau, Fp « 0.8. Thus, fluctuations play a significant role at all 
temperatures up to T = 1. However, above T « 0.4 the mean force is quite small and even the 
fluctuations cannot make it significantly large. As shown below, large force fiuctuations lead to 
significant fluctuations of the chain friction coefficient, which mark the failure of melt-based 
models. Well above Tq thermal noise overwhelms the activation energy, the forces become 
negligibly small and the dynamics should converge to the Tube model. 

To link these results to rheological behaviour, it is prudent to first establish the relation 
between chain friction and the disentanglement force distribution. Pulling a chain at a rate 
V consists of alternating uncoiling and disentanglement processes, giving rise to an oscillating 
pulling force. As the (n — l)th EP dientangles the segment following it uncoils behaving 
as an entropic spring whose effective length is the chain segment between the pulling point, 
n = 0, and the nth EP. Once the force reaches the disentanglement threshold of the nth EP 
it disentangles and the force drops. Presuming equilibration on the time scale of this process, 
the drop of the force is to the value consistent with elasticity of an entropic spring longer by 
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Fig. 4 - The relative width, F — {5f'^)^^'^/{f) as a function of temperature. 



the segment between the nth and (n + l)th EP. Now, the number of EPs disentangled per 
unit time r is iV = (vt/Io), where is the typical distance along the chain (measured in units 
of Kuhn's length a) between two successive EPs [|| . The force invested is 

F = r-^y, fn{t)dt = {f)M (7) 

where tn is the time between the disentanglements of the (n — l)th and nth EPs, fn{t) is the 
force during the time interval t„, and (...)jv indicates a time average for a given number of 
EPs N . To maintain a constant pulling rate, v, this force is balanced by friction and therefore 
= F. Therefore, fig. | gives (up to a numerical factor) ^(T) for a particular value of 
N . If the terms on the right hand side of cq. ^ have a typical value, g, then the friction 
is g/lo- If, however, /„ is broadly distributed then Eq. (Q) makes it immediately apparent 
that ^ can also suffer large fluctuations and its measured value may depend nontrivially on 
N . Furthermore, the sensitivity of the distribution of /„ to temperature affects directly the 
temperature dependence of the mean and the fluctuations of ^. We therefore expect ^ to be 
characterised well by a typical value only above Tq, where the distribution of /„ is narrow. 
Below To its measured value should appear noisy. Large fluctuations in ^ undermine melt- 
based models of polymer dynamics. For example. Rouse dynamics are based on a Langevin- 
type equation where the chain friction coefficient is presumed to have a typical constant 
value, and this description breaks down when the fluctuations of ^ become relevant. Similarly, 
the Tube model needs to be reconsidered when the chain friction is widely distributed. It is 
therefore suggested here that, in the absence of other mechanisms for the dynamics, Tg ~ 0.4 is 
the crossover temperature below which melt-based models for polymer dynamics break down. 
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Since one expects such models to fail around the glass transition then Tg « Tq, suggesting 
that Tg may be found from molecular data. For example, for polystyrene Tg « y!QK and 
therefore Eb/ks ~ 95QK and /o « 50pN. Thus this conjecture can be tested experimentally 
using the recent techniques of single chains pullout. 

To conclude, I have derived the distribution of forces required to disentangle non-permanent 
entanglement points as a function of temperature. I have shown that this distribution is nar- 
rowly peaked around small forces for temperatures in the melt and it broadens with cooling. At 
low temperatures a gap develops which corresponds to the increasing difficulty to disentangle 
at low pulling forces, as one would expect. The distribution, and consequently the dynamics, 
have been shown to depend on only two molecular quantities: Eh, the barrier energy needed 
to slide one chain across another, and e = Et,/ {a /scission)- The pulling force fluctuations have 
been found to dominate the dynamics below Tq ~ OAEb/ks and the implications of these 
results to viscosity on the molecular level have been discussed. This analysis pinpoint the 
range of validity of melt-based models by quantifyings how and at what temperature they 
break down. It also forms a basis to construct a consistent model that can bridge between 
the dynamics in the melt and glass states. The present description holds as long as the chain 
consists of both coils and non-pcrmnanent entanglement points and is not valid for too slow 
an equilibration or too fast a pulling rate. 

Since the chain dynamics are fundamental to most aspects of polymer systems these results 
are expected to have a wide range of applications. For example, they enable to compute 
directly the mean force required to pull an ensemble of chains along an interface in the interior 
of a polymer system, which is directly relevant to yield stresses of glassy polymers. Carrying a 
similar analysis on interfaces between different polymer materials (joined by welding, grafting, 
or any other method) would give the strength of such joints. An application to the relaxation 
behaviour of rubbers is being currently considered by this author. On a fundamental level, it 
is possibloe to re-analyse the Rouse dynamics with a noisy friction coefficient using the results 
obtained here, extending the usefulness of the traditional approach to lower temperatures 
than it is currently valid for. This analysis also paves the way to first principles calculations 
of the viscosity around the glass transition directly from molecular information. 

* * * 

It is a pleasure to acknowledge helpful discussions with Prof. Sir S. F. Edwards, Dr E. 
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